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Abstract 

We study a boundary layer problem for the Navier-Stokes-alpha 
model obtaining a generalization of the Prandtl equations conjectured 
to represent the averaged flow in a turbulent boundary layer. We 
solve the equations for the semi-infinite plate, both theoretically and 
numerically. The latter solutions agree with some experimental data 
in the turbulent boundary layer. 

1 Introduction 

Boundary-layer theory, first introduced by Ludwig Prandtl in 1904, is now 
fundamental to many applications of fluid mechanics, especially aerodynam- 
ics. 

Consider the case of two-dimensional steady incompressible viscous flow 
near a flat surface. Let x be the coordinate along the horizontal surface, 
y be the coordinate normal to the surface, and (n, v) be the corresponding 
velocity of the flow. Near the wall u is significantly larger than v. Also, 
u changes in y much faster than it does in the x coordinate. Then, in the 
case of a zero pressure gradient, the following Prandtl equations are used to 
approximate the Navier-Stokes equations. 

dx ^ dy W (I) 

where is a kinematic viscosity; density is chosen to be identically one. 
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A classical example of the boundary layer is a flow near the semi-infinite 
plate {{x,y) : x > 0,y = 0}. In 1908, Blasius discovered that in this case 
there exists a similarity variable C = Equations ([l|) can be reduced to 
an ordinary differential equation 

h'" + hh" = (2) 

with h{0) = h'{0) = 0, and h'{^) ^ 1 as ^ ^ oo. Then 

u{x,y) = Ueh' (^=] , v{x,y) = -^h' (^=] 

are solutions of (||) and they match experimental data in the laminar bound- 
ary layer. Here Ue is the horizontal velocity component of the external flow, 
k = —, and Rx = ^. 

For high local Reynolds numbers Rx the flow becomes turbulent and 
analogs of the equations (^) and (^) were not known. In this paper we 
use Navier-Stokes-alpha model of fluid turbulence also known as viscous 
Camassa-Holm equations to study averaged flow, and obtain turbulent bound- 
ary layer equations that generalize (||) and (^) . Our reduction to an ordinary 
differential equation (see (^)) also uses Blasius's similarity variable. 

Alpha-model was first introduced in Q as a generalization to n dimen- 
sions of the one-dimensional Camassa-Holm equation that describes shallow 
water waves. Later, NS-a model was proposed as a closure approximation 
for the Reynolds equations, and its solutions were compared with empirical 
data for turbulent flows in channels and pipes 

2 Derivation 

We study two-dimensional steady incompressible viscous flow near a flat 
surface. Let x be the coordinate along the surface, y be the coordinate 
normal to the surface. Denote also (n, v) to be the velocity of the flow. 
2-D Navier-Stokes-a model is used to study a boundary layer flow. 

J + (u • V)v + VjVuj = I'Av - Vq , . 

\ V • u = 0, ^"^^ 

where u = {u,v), and v 
way: 

V 



(7, r) is a momentum defined in the following 



" dxi \ dxj 
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Non-slip boundary conditions are used: u|y=o = 0. The other boundary 
conditions are going to be determined later. 

Fix I on the x-axis and define e(/) to be e := l/^/Ri = \J v j (uj,)-, where 
V is viscosity and Ug is the horizontal velocity component of the external 
flow. 

We change variables xi = \x, yi = ^y, m = ^n, vi = -^v, qi = ^g. 
ai = ^. In addition, assume that ai doesn't depend on y variable. After 
neglecting terms with high powers of e, dropping subscripts and denoting 

we derive the following turbulent boundary layer equations that generalize 
Prandtl equations: 

^^|^ = -*9 (4) 

Note that derivative of q is not zero in y-direction. Therefore, we introduce 
Q{x,y) := q + \v? — (^^) ' ^'^'^ obtain the following equations: 

|Q = o 
l^ + |^ = o- 

Consider now a case of a two-dimensional steady incompressible viscous 
flow near the semi-infinite plate {{x,y) : x > 0, y = 0}. The following 
assumptions are made in this case: 

• a = ^/x^i 

• Zero pressure gradient, i.e. = 0. 

In addition, we will study the solutions (tioo 5 ^00) of (0), which on some 
adequate interval xi < x < X2 are of the form 

1 y 

Jx \ X 



Let /i(^) = Jq f{ri)d'q. Then g = i^/i' — ^h, and we have the following 
equation for h: 

h'" + ^hh" - 0" (h'"" + ^hh"A = 0. (5) 
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The boundary conditions are h{0) = h'{0) = 0, h"{0) > 0, and h'{^) 1 
as ^ ^ oo. Note that if h{^) is a solution of (|5|), then h{x) := f3h{Px) is a 
solution of 

h'" + ^hh" - (^h'"" + Ihh""^ = 0. (6) 

This equation can be also written as m"' + ^hm" = 0, where m = h — h" . 
The following theorem is proven for this equation in Q. 

Theorem 2.1 Given any a > 0, b, there exists c{a,b) such that h'{^) — > 
const. > as ^ oo, where h is a solution of (d; with h{0) = h'{0) = 0, 
h"{0) = a, h"'{0) = b, h""{0) = c. 



3 Comparison with Experimental Data 

It is common to use 



= , = — 



in the turbulent boundary layer, where Ur = y ^ff 
zontal axis and denote = Rx = f^- As argued in the section |2|, 



. Fix X on the hori- 

y=0 



U = Uph! 



y 



represents a horizontal component of the averaged velocity for some h sat- 
isfying (|) with h{0) = h'{0) = 0, h"{0) = a > 0, h"'{0) = b, and h'{0 1 
as ^ ^ oo. 

For any such h, h{^) = l3h{f3^) satisfies (|). Then p"^ = \im^^^h'{^). 
In addition, 

Given c/, a skin- friction coefficient and Re, a Reynolds number based on 
momentum thickness we find a, b, and R^ so that the following conditions 
hold: 

1. Cf = 2/{miyU+f . 



2- Re = ll^u{l-if;)dy. 



3. Von Karman log law for the middle inflection point in logarithmic 
coordinates. 
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A family of curves {u}cj,_Rg was compared with experimental data of 
Rolls-Royce applied science laboratory, ERCOFTAC t3b test case (see Fig. 
|l|, |2|, and |3|). Comparison shows that the case a + (3b > corresponds to a 
laminar region of a boundary layer, a + fib < corresponds to a turbulent 
region of a boundary layer. The case a + (3b = corresponds to a transition 
point. 
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Figure 3: c/ = 4.01, Rq = 1436, /3 = 19.8 



